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r***- \ In this work, a model for the np — > (i(7r7r) reaction is developed. It is shown that 

the structure of the deuteron momentum spectra for a neutron beam momentum 
of 1.46 GeV can be explained as a consequence of the interplay of two mechanisms 
involving the excitation of the iV*(1440) resonance and its subsequent decay into 
^ ■ ^ ( 7r7T )'s-wave arm respectively. The relevance of the present analysis for the 

study of the Roper excitation and decay properties, as well as for the interpretation 
of other two-pion production experiments is discussed. 
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Abstract 



1 Introduction 



The study of nucleon-nucleon inelastic collisions provides a powerful tool to 
deepen our insight into the properties of nucleon-nucleon interactions and 
^ ■ baryonic resonances. A large amount of theoretical work on threshold meson 

production has been performed over the last years [1] , stimulated by the precise 
data obtained at IUCF, CELSIUS and COSY [2]. On the other side, double 
pion production reactions in the nucleon, 7 A — > Nirir and ttN — > Ntttc, have 
proved to be essential as a test of Chiral Perturbation Theory [3] and as a 
source of information about A* (1440) and A* (1520) [4-7]. In this context, 
the still scantily explored two pion production channel in the collisions of 
nucleons and light nuclei appears as promising research area. 

Most of the work, both experimental and theoretical, on two pion production 
in nucleon-nucleon collisions was performed in the seventies and in connection 
with the ABC effect. The ABC anomaly is an enhancement in the missing mass 
spectra close to the tttt production threshold, observed for the reactions pd — > 
3 HeX [8], np — > dX [9] and dd 4 HeX [10]. Although any interpretation of 
the ABC as a resonance is excluded [9], the origin of it is still poorly understood 
[9,11]. An important step towards the understanding of the ABC effect has 
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been taken in Ref. [12], where the 4 He spectra from the dd — > 4 HeX reaction 
at a deuteron beam energy of 1250 MeV [10] has been explained assuming 
that pions are independently produced in reactions involving two different 
pairs of nucleons from the projectile and target deuterons. Nowadays, two 
pion production in pp collisions is being studied experimentally at CELSIUS 
for beam kinetic energies between 650 and 775 MeV. On the theoretical side, a 
microscopic model for the NN NNnir reaction has been recently developed 
[13]. It includes mechanisms with the excitation of N* and A resonances, as 
well as some non-resonant contributions, and gives a satisfactory description of 
the available experimental data on total cross sections for most of the channels 
and in a wide range of energies up to 800 MeV above threshold. 

In this letter, I focus the attention on the deuteron spectrum in np — > d{itit) 
measured by Hollas and collaborators [14] using a nearly monokinetic neutron 
beam with central momentum p n = 1.463 GeV. This experiment is somewhat 
similar to the one of Plouin et al. [9] , but at lower energies ( T n = 795 MeV 
in [14] vs 1160 MeV in [9] ). Therefore, the analysis is simpler since, once 
the 7T° peaks are subtracted, only the double pion production mechanism is 
present. Apart from that, one expects that, being closer to threshold, the 
reaction mechanism might be simpler. The ABC peaks are not present in the 
data; they rather show a well defined bump at high mr missing masses, in 
disagreement with the models available in the literature [15,16]. From this 
comparison, the authors concluded [14] that neither double-A formation nor 
double-nucleon exchange provides the appropriate description of the reaction 
at p n = 1.46 GeV. A similar enhancement has also been observed for the 
reaction pd — > 3 Heir + 7r~, which is being studied using a beam of protons 
from COSY ( MOMO experiment ) [17], at a Q value close to the one of the 
experiment by Hollas et al. [14]. These common features point to a common 
dynamical description. 

Here, it is shown that the deuteron spectra for np — > d(7m)° at p n = 1.46 GeV 
can be understood as a consequence of the interference of two mechanisms 
involving the excitation of the Roper resonance iV*(1440) and its subsequent 
decay into N (^ir)'s-Zave an d A7r respectively. The implications of this finding 
to the MOMO experiment are also discussed. 



2 The model 

The model is schematically presented in Fig. 1. It is a reduced version of the 
model of Ref. [13], modified for the case where one has a deuteron instead of 
two free nucleons in the final state. The choice of the mechanisms was based 
on their contribution to the total cross section for the pn — > pnix+ix~ reaction; 
the situation for the pn — > pnir°-K channel is similar. At T p = 800 MeV, 
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the mechanism with N* — ■> N^ittt)^^ ave gives cr ~ 11 fib, being by far the 
most important. The second largest contribution comes from iV* — > A7T with 
cr ~ 0.5 fib; as we will see, in the case of a deuteron in the final state, its 
contribution is larger with respect to the dominant iV* — > iV (^)^Z^ ave and 
crucial to obtain the right shape. All other mechanisms give a < 0.3 fib; I do 
not include them all, but just the double- A ( a ~ 0.1 fib ) one in order to 
make contact with the model of Ref. [16]. 

The deuteron momentum spectrum is the sum of the 7r + 7r~ and the 7r°7r° 
contributions. In the Laboratory frame, the charged pions piece is given by 



d 2 a 1 1 MM d {p' d ) 



dp' d dn' d 4(2tt)5 E> d p nPd 



J dE n dip n - \M Rri r 2 \ 2 \ ■ (1) 

V Rrir 2 / CM 



Here, E' d and p' d are the deuteron energy and the modulus of its momentum, 
both in Lab. frame; p d is the modulus of the deuteron momentum in the 
center-of-mass system (CM); M and M d stand for the nucleon and deuteron 
masses respectively. The integral in brackets must be calculated with all the 
kinematical variables defined in CM; it runs over the polar angle and the 
energy of one of the outgoing pions. The amplitude squared is summed over 
the deuteron spin (R) and averaged over the spins of the incoming nucleons 
(ri, r 2 ). For the neutral pions channel, the expression is the same but with 
and extra 1/2 factor, which is a consequence of having two identical pions in 
the final state. The difference of masses between charged and neutral pions is 
taken into account for the phase space, but isospin symmetry is assumed in 
the calculation of the amplitude. 

The amplitude can be expressed as 

M Rrir2 = Q r ^ r 2 
{((pr[\ Vi \pn) (nr' 2 \ V 2 \nr 2 ) - (nr[\ Vy \pn) (pr' 2 \ V 2 |nr 2 )) tp d (P 2 ) 

+ {{P r l\ \Wl) ( nr 2\ Vl \ nr 2) - { nr l\ V2 \Wl) (P r 2\ Vl l nr 2)) ^ci(Pl)} (2) 

where y3d(P) is the s-wave deuteron wave function in momentum space, nor- 
malized as 

/t^HI. (3) 
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dq 
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T=0,1 
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The d-wave part has been neglected. For the wave function, different expres- 
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sions and parameterizations can be used [18-20]. The value of Pi (2) depends 
on the mechanism; for those of Figs, la and lb, Pi (2) — q + Pi(2) — Pd/2 and 
for the A - A mechanism ( Fig. lc ) Pi( 2 ) = q + Pi(2) - Pd/2 - p^, Pi(2) 
and Pd been the momenta of the proton (neutron) and deuteron respectively; 
p,,- is the momentum of the pion, over whose energy the integral in Eq. 1 is 
performed. The function Dt=o,i{q), which stands for the meson propagators 
and form factors, will be discussed later; q = (go, q) is the four momentum 
transfer from one nucleon to the other; q is given by energy conservation 
in the vertices and, therefore, depends on the energy of one of the outgoing 
nucleons, taken to be one half of the deuteron energy. 

The matrix elements in Eq. 2 are evaluated for the different mechanisms using 
the Feynman rules that can be obtained using phenomenological Lagrangians; 
some of the required ones are 

Un, = —^ A S}(d i( f>)T^ N + h.c. , (4) 

£n*Nit = — iplf*<Ti(di(t>)Tip N + h.c. , (5) 
= g -^As\{di<t>)T^ N * + h.c. . (6) 



In Eqs. 4, 6, (T^) are the spin (isospin) 1/2 — > 3/2 transition operators [21]; 
ipNi V'a, ipN* and <f> stand for the nucleon, Delta, Roper and pion fields, while 
m n is the pion mass. The absolute value of the coupling constants /* = 2.13, 
/ = 0.477 and Qn'A-k = 2.07 are obtained from the partial decay widths of 
the A and iV*(1440) [22]. In the case of the decays N* -> Ntt and N* -> Avr, 
branching ratios of 65% and 25% respectively are assumed, as well as an iV* 
total width of 350 MeV; the signs correspond to those provided in earlier 
analyses of the (n, nn) reactions [5,6]. The Lagrangian for the iV7V7r vertex is 
the standard one which, in the non-relativistic limit, looks like 

£nntt = ^^ L il>N<7i(d i <l>)Trl) N , (7) 

with f NNw = 1. 

A general Lagrangian for the N* — > N (tctt)^^^ decay [23] is 
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(d^d v i) N *)(Td^(j))(Td u (t))i) N + h.c. 



(8) 
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where / = 92.4 MeV is the pion decay constant and M*, the mass of the 
Roper resonance. Using the partial decay width, the parameters c* and c* 2 
can be constrained to an ellipse [23]. In order to further constrain them, the 
model of Ref. [6] and the overall data on 7r~p — > 7r + 7r~n have been used [24]. 
Assuming a branching ratio of 7.5% for the N* — > N (^)'sZ t ^ ave decay, the 
best agreement is obtained for c* = —7.3 GeV^ 1 and c* 2 = (Set I), but the 
data seem to be still compatible with the choice of c* = —12.7 GeV -1 and 
(? 2 = 2.0 GeV -1 (Set II). In this study, Set I will be used except where a 
different choice is explicitly indicated. 

The matrix elements of Eq. 2 contain, apart from the vertices described above, 
Roper and Delta propagators, given by 

A(p) = . 5 * i T , w = ' = ; l = (N*,A) (9) 

with p = (po,p) the momentum of the resonance and Ti(p), its total width. 
The partial decay A — > Nn practically accounts for the total Delta width. In 
the case of the Roper, the major part of the width comes from the decay into 
nucleon and pion and the rest from the decay into nucleon and two pions. All 
of them, except the small (less than 8% at the N* peak) N* — > Np decay, can 
be calculated with the Lagrangians given above. 

Finally, let us consider Dr=o,i(<?)- For the T = 1 potential, I calculate the 
diagrams assuming a pion exchange and make the substitution 

DP =1 {q)m - VL(Q)QiQj + ^4(<?)(<% - W (10) 

so that Dt=i includes the longitudinal pion exchange, the transversal rho 
exchange, and the short range correlations that take into account the repulsive 
force at short distances. Functions V[(q) and V^(q) are described elsewhere 
[13,25]. With respect to the T = channel, in a recent analysis [26] of the 
(a, a') reaction on a proton target, the strength of the isoscalar NN — > iVW* 
transition was extracted by parameterizing the transition amplitude in terms 
of an effective "a" which couples to iViV as the Bonn model a (s^aw/^ 71 " — 
5.69) [19] and couples to iViV* with an unknown strength provided by a best 
fit to the data {g 2 rNN */An = 1.33). According to the fit, g a NN and g a NN* have 
the same sign. Therefore, we have 

1 /A 2 -m 2 \ 2 
q l — m~ \ Aj — q 1 J 

where an equal A CT = 1.7 GeV for the form factor is assumed for both vertices. 
The explicit expressions for the amplitudes are too involved to be reproduced 
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here, but some details of their evaluation can be traced in Ref. [13], where all 
the amplitudes for the free reaction pp — > ppn + n~ are given in the Appendix. 



3 Results and discussion 

With the ingredients described in the previous section, one can calculate the 
deuteron spectra at p n = 1.46 GeV for different angles. The results are shown 
in Fig. 2. They compare quite well with the data of Hollas et al. [14], and 
certainly much better than the previous models [15,16]. The curves, in general, 
underestimate the data, maybe because there are many other mechanisms, 
not considered for the sake of simplicity, that, even been individually small, 
could in sum enhance the cross section. Some of the approximations made, 
like the neglect of the deuteron d-wave and non-relativistic approximation in 
the vertices, as well as the intrinsic uncertainties of the mechanisms included 
can also be a source of discrepancies; they are discussed below. The large data 
point at the edges of the spectra show the contamination of the ir° peaks 
[14]. I have also estimated the influence of the width of the neutron beam 
by averaging the double differential cross sections over a Breit-Wigner profile 
of 40 MeV width and centered at p n = 1.46 GeV [14]; the contribution of 
this effect is very small and irrelevant for our study. The plotted curves are 
obtained using the deuteron wave function derived from the Paris potential 
[18]; with the Bonn wave functions [19] the results are very similar, while for 
the Hulthen one [20] the distributions are overall larger, up to a factor two at 
the position of the central bump. 

The mechanism N* — > iV (^)'s-wave (Fig- l a ) produces spectra very similar 
to phase space; its contribution is certainly the largest, but its relative size 
with respect to the N* — > An mechanism is not as large as one would naively 
expect from estimations based on the total cross sections obtained for both 
mechanisms in the free NN — > NNnn reaction. Nevertheless, it is not sur- 
prising since here we are sensitive only to a reduced phase space region and 
choosing a particular combination of the quantum numbers of the outgoing 
nucleons (those of the deuteron). For this mechanism, I have calculated the 
contribution of the d-wave part of the deuteron wave function and found it 
negligible. 

The N* — > An mechanism (Fig. lb) exhibits a wide bump at the center of 
the spectra (high nn masses), whose maximum falls fast with the increase 
of the deuteron angle, and small peaks at the edges of the spectra (low nn 
masses); the size of these peaks does not vary appreciably with the angle. This 
mechanism plays a crucial role in providing the right shape to the distributions 
through its interference with the larger NN — > NNnn contribution. This 
interference is constructive at high nn masses and destructive at low ones. 
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Such pattern can be understood by realizing that the N* — > An amplitude is 
dominated by terms proportional to the scalar product of the outgoing pions 
three momenta; this scalar product has different signs in the center of the 
spectra, where the pions go back to back, and at the edges, where they travel 
together. In order to further illustrate the effect of the interference, Fig. 3 
shows the effect of changing the relative sign of the two amplitudes. The data 
clearly favor a choice of the sign of gN*An in agreement with earlier works [5,6]. 

The double- A mechanism (Fig. 1) is so small that it can hardly be distin- 
guished in Fig. 2. In Fig. 4, the contribution of this mechanism alone is shown 
for p n = 1.5 GeV, Q\ a b = 4.5° and using the Hulthen wave function, in order 
to compare with the result of Bar-Nir, Risser and Shuster (Fig. 4 b of Ref. 
[16]). The differential cross section obtained in the case of only pion exchange 
is very similar to the one given by the relativistic model of Ref. [16]; the inclu- 
sion of the rho exchange modifies the result, and the short range correlations 
between the initial nucleons cause a strong reduction of the strength of this 
mechanism. At T p = 775 MeV (p p = 1.43 GeV), the pp — > pnnn and pp — > dnn 
reactions are probably dominated by the double-A mechanism; therefore, the 
future data from CELSIUS would provide an excellent opportunity to study 
the interplay of the described ingredients and, in particular, the role of the 
short range correlations. 

Wherever dealing with the Roper resonance, the lack of a precise determina- 
tion of its properties is a problem that should be faced. In Fig. 5, I investigate 
how the results vary with some of the uncertainties. As can be noticed, the 
shape is not affected by these changes, but some set of values are preferred. 
Fig. 5 a, shows how the spectrum changes with the modification of c\, c% within 
the previously accepted values, while keeping the total width and the partial 
branching ratios fixed; set II gives a better agreement with the data. In Fig 5 
b, we show the range of uncertainties that come from the variation of the total 
width of the N* in the limits given by the Particle Data Book [22], that is from 
250 MeV (lower line) to 450 MeV (upper line). Finally, the dependence on the 
partial branchings of the N* to N^n)^^^ and An, with the total width 
fixed to 350 MeV, is studied in Fig. 5 c. The dashed line corresponds to a 5% 
branching of N* to N(nn)^z^ ave and a 30% to Air, while for the dash-dotted 
one, a larger 10% of N(nn)^w aV e an d a smaller 20% of An; the data prefer 
the latter choice. These data alone do not allow to disentangle the different 
effects, but, in combination with other data that will be available in the future 
(like the data on pp — > dnn at T p = 600 — 775 MeV from CELSIUS), would 
be an important source of information about the Roper resonance. 

As mentioned in the Introduction, the present model could provide an expla- 
nation to the deviation from phase space observed in the two pion invariant 
mass distribution of the reaction pd — > 3 Hen + n~ , studied with a proton beam 
of momentum p p = 1.15 GeV [17]. This experimental result has been inter- 
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preted with the hypothesis that the pions are preferably emitted in p-wave 
[27]. An ansatz based on this idea can explain the data reasonably well but, 
unfortunately, the initial tail of the p meson produced via pd — * 3 Hep can 
hardly produce the required strength. In Fig. 6, we show the same observable 
for the reaction np — > dnir; the shape obtained for the charged pions chan- 
nel is very similar to the one observed at COSY. The neutral pion channel, 
though, exhibits a shape similar to the one of the charged pions, and only a 
factor about 1/2 smaller at the peak position, in apparent contradiction with 
the much smaller rate of n° pairs, compared to the charged ones, obtained in 
the MOMO experiment. This is a consequence of the fact that, in the present 
model, the pions are almost always produced in s-wave. However, we should 
bear in mind that the interference pattern can be very different when the 
nucleons are bound in a 3 He nucleus instead of a deuteron. Detailed calcu- 
lations are, therefore, required to check if the model can explain the MOMO 
experiment. 



4 Conclusions 



In summary, a simple model for the np — > d(ii7r) reaction has been developed, 
based on a previous model for the free NN — > NNnir and including the most 
important resonance contributions. It is shown that the bump in the center 
of the deuteron momentum spectra (high nn masses) observed at a neutron 
beam momentum of p n = 1.463 GeV (T p = 795 MeV) can be explained as 
a result of the interference of two mechanisms involving the excitation of 
the Roper resonance: the dominant and phase-space like N* — > N (irir)^^^ 
(Fig. 1 a) and the smaller in size N* — > An (Fig. 1 b), but determinant to 
obtain the right profile. The mechanism of double-A (Fig. 1 c) excitation, 
considered in an earlier model for the same reaction, but only including pion 
exchange, is significantly reduced by the short range correlations. Other two 
pion production reactions like pp — > prmix and pp — > dmr, currently studied 
experimentally at CELSIUS for energies between 650 and 775 MeV, would be 
helpful to clarify this issue, The size of the spectra depends appreciably on the 
Roper resonance and can be used, together with other two pion production 
reactions, to learn more about this resonance and its decay properties. Finally, 
the present model could be helpful to understand the two-pion invariant mass 
distributions observed from pd — > 3 He7r + n~ at COSY. 
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Fig. 1. Set of diagrams of the model. 
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Fig. 2. Deuteron momentum spectra for np — ► d(-7r7r)° at p n = 1.46 GeV and different 
laboratory angles (solid lines) compared to the measured data [14]. The dotted line 
corresponds to the N* — > N(7nr)sZS;at;e mechanism (Fig. 1 a); the short-dashed line 
stands for the N* — > A-zr (Fig. 1 b) and the long-dashed one, for the double-A (Fig. 
lc). 
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Fig. 3. Calculated spectra for two different choices of the gN*An sign. The data 
clearly favor the positive sign (solid line) with respect to the negative (dash-dotted 
line) . 



p„ = 1.5 GeV 
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Fig. 4. Contribution of the double-A mechanism in the case of tt exchange alone 
(dashed line), ir + p exchange (dash-dotted line) and tt + p+ short range correla- 
tions (solid line). In this case the calculation uses a Hulthen wave function for the 
deuteron. 
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Fig. 5. Dependence of the spectrum on certain decay properties of the TV* (1440). 
The solid line in (a) shows the result for the Set I of parameters c*, c* 2 used ev- 
erywhere else (c| = 0), and the dash-dotted one is obtained for Set II. In Fig. (b), 
different values of the N* total width are considered: dashed line, 250 MeV; solid 
line, 350 MeV; dash-dotted line, 450 MeV. The branching ratios of the different de- 
cay channels of the ./V* are modified in (c), with the total width fixed to 350 MeV; 
dashed line: ^(A^Tr)^^) = 5% and Sr(Avr) = 30%; solid line: 7.5% and 25%; 
dash-dotted line: 10% and 20%. 
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Fig. 6. Two-pion invariant mass spectrum for np — ► dirir at p n = 1460 MeV for both 
charged and neutral pions. 
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